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ABSTRACT 

We present a spectroscopic sample of 910 distant halo stars from the Hypervelocity Star survey from 
which we derive the velocity dispersion profile of the Milky Way halo. The sample is a mix of 74% 
evolved horizontal branch stars and 26% blue stragglers. We estimate distances to the stars using 
observed colors, metallicities, and stellar evolution tracks. Our sample contains twice as many objects 
with i? > 50 kpc as previous surveys. We compute the velocity dispersion profile in two ways: with a 
parametric method based on a Milky Way potential model, and with a non-parametric method based 
on the caustic technique originally developed to measure galaxy cluster mass profiles. The resulting 
velocity dispersion profiles are remarkably consistent with those found by two independent surveys 
based on other stellar populations: the Milky Way halo exhibits a mean decline in radial velocity 
dispersion of —0.38 ± 0.12 km s^^ kpc^^ over 15 < i? < 75 kpc. This measurement is a useful basis 
for calculating the total mass and mass distribution of the Milky Way halo. 
Subject headings: Galaxy: halo — Galaxy: kinematics and dynamics — 



1. INTRODUGTION 

A fundamental property of the Milky Way galaxy is 
its total mass. In the Cold Dark Matter paradigm, the 
total mass of a galaxy's dark matter halo correlates with 
its merger history, star formation history, and number 
of its satellite sub-halos. A galaxy's mass distribution 
is also a fundamental constraint on theories. Cold dark 
matter model s predict that densit y follows an universal 
NFW profile (jNavarro et al.lll997f l: in Modified Newto- 
nian Dynamics there is no dark matter halo and the mass 
distribution is highly flattened compa red to cold dark 
matter models fjHernandez et al.ll2009 i). The Milky Way 
provides an unique laboratory for testing these issues. 

Nevertheless, the mass and mass distribution of the 
Milky Way are among the most poorly known of all 
Galactic parameters, especially at large distances i? > 50 
kpc. Total mass estimates span at least a factor of 4, 
from 0.5 x 10^^ to 2 x 10^^ jvIq (see below). Recent 
theoretical work based on semi-analytic models and tim- 
ing argur nents suggest the Milky Way m ay be even more 
massive (iLi fc White! [2008t lAbadi ellll [20091. Clearly, 
the mass and mass distribution of the Milky Way remain 
controversial and important issues. 

A powerful approach for measuring the Milky Way's 
mass distribution is measuring the motions of tracer ob- 
jects. Radio m asers probe the G alactic rotation curve 
out to 15 kpc (jReid et al.l l2009( l: H l g as probes the 
rotat ion curve to larger distances (|Kalb erla & Dcdes 
|2008( ). Luminous tracers at i? > 50 kpc, however, 
are rare. Historically, globular clusters and dwarf 
galaxies were u sed to measure the total mass of the 
Milky Way (e.g. 'Little & Tremaine"1987'; 'Zaritsky et alj 
1989; Kulessa &LLvnden-BcU 1992; Kochanck 199Q; 



Wilkinson fc Evansi Il999l ). These mass estimates were 
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based on samples of a few dozen objects, and suffered 
from a systematic factor of ~2 uncertainty depending 
on the inclusion of Leo I. Post-main sequence halo stars 
provide denser tracers , but the blue horizontal br anch 
(BHB ) star sain p les of lSommer-Larsen et al.) ()1997t ) and 
Saka moto et alJ (l2003f) are limited to i? < 20 kpc. 
iBattaglia et al.l (|2005l ) added 58 distant red giants from 
the Spaghetti survey and claimed a large dec hne in the 
velocity dispersion at R >50 kpc. Recently, IXue et ahl 
(jl308)' analyzed the Sloan Digital Sky Survey (SDSS) 
sample of 2,466 BHB stars and found a small decline in 
velocity dispersion with distance. Only 80 of the SDSS 
BHB stars are located at i? > 50 kpc. 

Here we use the distant halo stars from the hyperve- 
locity staxJJTVS) program (Brown ct al. 2005a, 2006a, b|, 
I2007bll3 . [2009al fll) to measure the velocity dispersion pro- 
file of the Milky Way. Our dataset is a complete spectro- 
scopic sample of 910 stars observed over 7300 deg^ of the 
SDSS Data Release 6 imaging region. We gain increased 
leverage on the velocity dispersion profile for i? > 50 kpc. 
We find an average decline in radial velocity dispersion 
of 0.38 ± 0.12 km s'^ kpc'^ over 15 < i? < 75 kpc. 

In §2 we describe the observations and luminosity esti- 
mates for stars in our sample. In §3 we make parametric 
and non-parametric estimates of the velocity dispersion 
profile. We also discuss possible systematics, and com- 
pare our results with earlier work. We conclude in §4. 
The data are in the Appendix. 

2. DATA 

The hypervelocity star program is a radial velocity sur- 
vey of stars selected with the colors of late B-type stars. 
The radial velocity survey is now 93% complete over 7300 
deg^ of the SDSS DR6 imaging footprint. Here we focus 
exclusively on the stars with late-B and early- A spect ral 
types; we exclude all white dwarfs ()Kihc et al.ll2007airbl ). 
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Fig. 1. — Color-color diagram showing t he distribu t ion o f 
the HVS survey stars (solid squares) and the IXue etaLl 1I2OO8I) 
BHB stars (dots) compared to line s of constant absol ute magni- 
tude Mg for BHB stars ( red lines. IDotter et al.l 120081 ) and main 
sequence stars (blue lines IGirardi et al.|l2004|). All tracks are for 
[Fe/H]=— 1.9, the mean metallicity of halo stars. Dashed lines 
indicate the direction of increasing luminosity. At [Fe/H]=— 1.9, 
BHB stars and blue stragglers share identical luminosities around 
(u - g)o ~ 0.6. 



B s upergiants (iBrown et al.ll2007af). emissi on line galax- 
ies (iKewlev et al.ll2007l:lBrown et al ll2008bf) . and quasars 
(|Brownet al.ll2009aD . 

The dataset contains 910 stars: 571 stars from the orig- 
inal HVS survey (iBrown et al.iii2007G') . 331 stars from the 
new HVS survey ()Brown et al.ll2009aD. and 8 B HB stars 
from the earliest sample ( Brown et al.ll2005a^ . We be- 
gin by describing the observables - magnitude, position, 
velocity - all of which are well-determined. Stellar lumi- 
nosity is less well-determined and must be inferred from 
colors and metallicity. We discuss the luminosity esti- 
mates and distance determinations in some detail. The 
observed and derived quantities for each star are listed 
in Table lAll described in the Appendix. 

2.1. Photometry 

All photometry comes fr o m S PSS Data Release 
6 (lAdclman-McCarthiLeLin l2008f ). We use uber- 
calibrated PSF magnitudes, and c orrect the mag i iitude s 
and colors for reddening following ISchlegel et al.l (|1998f ) . 

2.2. Target Selection 

The HVS survey target selection emphasizes outliers 
in the halo population: we target stars redder in {u — 
g)o than known white dwarfs and bluer in (g — r)o than 
known BHB stars (Brown ct al. 2006b). This color cut 
through the stellar population, illustrated in Figure [H 
allows us to detect hypervelocity stars efficiently. The 
majority of targets, however, are normal halo stars. 

Our target selection includes stars with 1 7 < go < 19.5 
in the range —0.39 < {g — r)o < —0.25 (jBrown et al.l 
l2007c( ) and fainter stars with 19 < 50 < 20.5 over 
a broader color ra nge —0.40 < (5 — r)o < —0.20 
(jBrown et al.ll2009af) . We also include 8 confirmed BHB 




Fig. 2. — Distribution of velocities in the Galactocentric rest 
frame. The dotted line shows a fiducial Gaussian with zero mean 
and 106 km s~^ dispersion. The positive velocity outliers with 
Vrf > -1-400 km s~^ are the unbound HVSs. 



stars from IBrown et al.l (|2005af) with 19.5 < go < 20.25 
and —0.3 < {g — r)o < —0.1. All 910 targets are located 
in the SDSS DR6 footprint and have an average surface 
density on the sky of 0.12 deg~^. 

2.3. Radial Velocity 

We obtained spectroscopic observations at the 6.5m 
MMT telescope with the Blue Channel spectrograph. We 
operated the spectrograph with the 832 line mm~^ grat- 
ing in second order, providing wavelength coverage 3650 
A to 4500 A and a spectral resolution of 1.2 A. We ob- 
tained all observations at the parallactic angle, with a 
comparison lamp exposure for every survey object. 

We processed the data using IRAF^ in the standard 
way. We measure radial velocities by cross-corre lating 
the o bservations with radial velo city standards (I Fekel 
mil) using the package RVSAO (iKurtz fc Minklll99af) . 
The average radial velocity uncertainty of the stars is 
±12 km s-i. 

All velocities discussed here are in the Galactocentric 
rest frame, indicated Vrf- Given the outer halo location 
of the stars, we note that the observed radial velocities 
are almost purely (>85%) radial in the Galactocentric 
frame. We transform heliocentric velocities (vheiio) into 
Galactocentric rest frame velocities assuming a circular 
velocity of 220 km s^^ and a solar motion of (U , V, W) 
= (10, 5.2, 7.2) km s^^ (jPehnen fc Binnevlll998l ): 

Vrf = w?ie/io+220sinZcos5-|-(10cosZcos 6-1-5. 2sin^cos 6-1-7. 2sin6). 

iReid et all ([2009h argue for a larger circular velocity of 
250 km s~^ based on trigonometric parallaxes to star 
formation regions in the disk. We test using a circu- 
lar velocity of 250 km and find statistically identical 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 3. — Metallicity d istribution function of inner halo BHB 
stars IIBrown et al.ll2"008al) and outer halo BHB stars llXue et al.l 
120081) based on Ca II K. The shift t owards lower metallicity in 
the outer halo is expected; the IXue et al., (.2008.1 BHB stars with 
go > 17 have a mean [Fe/HJca = —1.9. 



velocity dispersion profiles. This insensitivity to circu- 
lar velocity arises because the stars in our high latitude 
survey have a mean value of |sin^cos6| = 0.32; chang- 
ing the Sun's circular velocity by 30 km s~^ results in 
a ±10 km s^^ change in the stars' rest frame velocities. 
10 km s^^ is smaller than our measurement error and an 
order of mag nitude smaller than the vel ocity dispersion 
of the star s iMcMillan fc BinnevI ()2009[ ) show that the 
iReid et all (|2009D data are consistent, at the 1-a level, 
with the canonical circular velocity of 220 km s~^. Thus 
we use 220 km s""'^ here. 

Figure[2]plots the resulting Galactic rest frame velocity 
distribution of the 910 halo stars. For reference, we also 
draw a Gaussian with zero mean and 106 km s~^ disper- 
sion (dotted line). The observations reveal a significant 
asymmetry of outliers in the tails of the velocity distribu- 
tion. There are no stars with Vrf < —300 km s^^ but 18 
stars with Vrj > +300 km s^^(thc HVSs). We address 
this issue in Section ??. 

2.4. Metallicity 

The strongest metal line in our spectra is the 3933 A 
Ca II K line. Unfortunately, at the effective tempera- 
tures sampled by the survey, 10, 000 < T^ ff < 15, 000 K 
the eq uivalent width of Ca ii K is small (i Wilhelm et al.l 
Il999al ). Thus, Ca ii K provides poor leverage on the 
metallicity of our stars. Metallicity is better determined 
for red der (cooler) BHB sta rs, for example the BHB sur- 
vey of iBrown et al.l (|2008af l and the BHB sample from 
the SDSS spectroscopic survey ()Xue et a l. 2008). The 
metallicity distributions of these two BHB samples are 
plotted in Figure [3] and are similar in shap e. How- 
ever, the outer halo sample of lXue et al.l ()2008f ) is about 
0. 2 dex more metal-po or than the inner halo sample 
of IBrown et"al] (|2008a^ . The mean metallicity of the 
IXue et al.l (|2008l ) BHB stars with 50 > 17 is [Fe/H]ca = 
-1.9. 

To make luminosity estimates, we assume our sur- 



vey stars have t he metallicity distribution function of 
IXue et aP ()2008f ). This assumption is reasonable given 
the very similar sky coverage and penetration into the 
halo of the two surveys: the stars occupy similar regions 
of the Milky Way halo. The reddest stars in our sam- 
ple, where we can e stimate metallicity, are metal-poor 
(|Brown et al.ll2006bf ) , consistent with the metallicity dis- 
tribution function of IXue et all (|2008D . 

2.5. Spectroscopic Identification 

Although the stars in our survey have the spec- 
tral types of late B- and early A-type stars, their 
nature is ambiguous. The old stellar population of 
the halo contains both evolved BHB stars and main 
sequence blue stragglers. Halo surveys consistently 
find that ^ -^50% of A-type s tars i n the field are blue 
stragglers (iNorr is & Hawki nj|1991l: [K inman e t al.lll994t 
' Preston et al."1994: Wilhclm et al.''l 999bc iClewlev et al l 




2004; Brown et al. 2003, 2005bl. l2008aHXue et al ' 



This result is problematic because BHB stars 
and blue stragglers can have very different lumi- 
nosities (Figure [T]). Spectroscopic measures of sur- 
face gravity can di scriminate the evolutionary s tate of 
the stars ( Kinma n et all 11994 IWilhelm etall Il999al : 
IClewlev et al.l 1200202001 . Unfortunately, surface grav- 
ity measures fail for our sample because BHB and main 
sequence stars have nearly identical surface gravities at 
the effective temperatures of our stars. 

We target stars so blue, however, that we largely ex- 
clude the possibility of blue stragglers. A ~0.75 Mq star 
with [Fe/ H]=-1.9 has a main sequence lifetime of a Hub- 
ble time (jCirardi et al.ll2004h . A field blue straggler can- 
not plausibly have more than twice the mass of a main- 
sequence turnoff star. In our sample, 36% of the stars are 
bluer than both (5 — r)o = —0.27 and {u — g )o = 0.87, the 
color of a 1.5 Mq star with [Fe/H]=-1.9 (jGirardi et all 
l2004f ). Thus the nature of many of the stars is clear: 
they are hot BHB stars. 

Stars bluer than a 1.5 M0 star are BHB; stars redder 
than a 1.5 Mq star are equally likely to be BH B stars 
or blu e stragglers. We base this conclusion on Xu e et al.l 
(j2008l ). who observe a BHB fraction of 47% at the red 
end of our sample. Overall, our sample is 74% BHB stars 
and 26% blue stragglers. 

2.6. Luminosity 

We estimate luminosity from stellar evolution tracks. 
iGirardi et all (12001120041 ) provide main sequence tracks 
in the SDSS passbands for metallicities ranging from so- 
lar to [ Fe/H]=— 2.3, but do not include the horizontal 
branch. iDotter et all (j2007, 2008) provide BHB tracks in 
the SDSS passbands over the same range of metallicities. 
To simplify the luminosity estimates, we fit low-order 
polynomials to the tracks as a function of both color and 
metallicity. 

We estimate luminosity two ways, using a star's {u—g)o 
and (g — r)o color. Our final luminosity is the weighted 
average of the two estimates. We weight by the slope of 
the observed {u — 17)0 versus {g ~ r)o distribution (Fig- 
ure [T]). This weighting is important because luminosity 
is most sensitive to (u — g)o at the blue end of our sam- 
ple and to [g — r)o at the red end of our sample. The 
weighting favors luminosities derived from (u — g)Q for 
(5 ~ ^)o < —0.24 and luminosities derived from (g — r)o 
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Fig. 4. — Distribution of luminosity Mg and velocity for fiducial 
stars in the bluest (left panel) and reddest (right panel) quartiles 
of our sample. The distributions are obtained by propagating the 
observed uncertainties of a go = 19 star at the median depth of 
our survey (au—a = 0.05, (Jg—r = 0.03, (t„ = 12 km s~^) with 
the lXue et al.l II2008I ) halo metallicity distribut ion function through 
stellar evolutionary tracks for main sequence l IGirardi et al.ll2004l) 
and BHB (Dottcr et al. 2003) stars. Blue stars have luminosity 
precise to 10%, whereas red stars have a bimodal distribution in 
luminosity (compare with Figure [TJ. 



for {g — r)o > —0.24. Figure [T] shows the resulting hnes 
of c onstant lumin o sity f or a star with [Fe/H]=— 1.9 from 
the lGirardi eTa P ('2004^ main sequence tracks (in blue) 
and from the [Dotter et al. (2008) BHB tracks (in red). 

The accuracy of the luminosity estimate depends on 
the accuracy of the stellar evolution models. For ex- 
ample, alpha-enhanced tracks (appropriate for the old 
stellar halo populations) are systematically bluer than 
solar-scale d tracks for both main sequence and BHB stars 
(|Lee et al.1 '2009). 

The precision of the luminosity estimate is more ro- 
bust because we apply the same tracks to all stars. The 
precision depends on observational uncertainties and is 
straightforward to quantify. A (70 = 19 star at the me- 
dian depth of our survey has uncertainties a(^u-g) = 0.05 
and a(^g_j.) = 0.03 in color and Oy = 12 km s~^ in ve- 
locity. We assume that m etallicity is randomly drawn 
from the 'Xue et aTl (|2008[ ) metallicity distribution func- 
tion (Figure [3|) . We then propagate these uncertainties 
through the main sequence and BHB tracks to visualize 
the resulting distribution of luminosity estimates. 

Figure [4] plots the distribution of luminosity and veloc- 
ity (centered at zero) for two fiducial stars in the bluest 
quartile (left panel) and the reddest quartile (right panel) 
of the survey. Stars in the bluest quartile have lumi- 
nosities precise to 10% because main sequence and BHB 
stars have essentially identical luminosities at these col- 
ors (Figure [1]). Stars in the reddest quartile, however, 
have bimodal luminosity distributions with a factor of 
4 spread in luminosity (a factor of 2 in distance). This 
degeneracy is broken for confirme d BHB stars: 33 o f our 
stars are listed as BHB stars in IXue et all (l2008f) and 
8 are BHB stars from our earliest sample ( jBrown et al.l 
[2005a). For these confirmed BHB stars we use only the 
BHB luminosity to calculate distance. 

Three of the stars with Vrf > -1-400 km s^^ in 
this surve y are confir r ned y oung main sequence 
B stars (^Fuentes et al.l I2006t iBonano s et al.' 120081: 
Lopez-Mor alcs & BonanosI 2008[ fPrzybill a et al.l 
2008allbl ). For the purposes of this paper, however, 
we treat these HVSs the same as the other stars in our 
sample: we estimate distances to the HVSs as if they 
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Fig. 5. — Density distrib ution of the HVS survey stars compared 
with the IXue et alj 1I2OOI) BHB sample. Our survey is complete, 
but we find an artificially shallow density profile because we observe 
over a broader color range at large depths go > 19 mag. 



were halo BHB stars or metal-poor blue stragglers. We 
clip the HVSs from the sample when calculating the 
velocity dispersion in Section ??. 

2.7. Distance 

We calculate distances from the observed apparent 
magnitudes and estimated luminosities: 



(2) 



where d is the heliocentric distance in kpc, go is the ap- 
parent magnitude corrected for extinction, and Alg is the 
absolute magnitude estimated above. We convert helio- 
centric distance d to Galactocentric distance R with the 
Sun at i? = 8 kpc. 

Figure [5] plots the resulting di stance dist r ibutio n of the 
survey stars compare d with thelXue et aD (|2008f l sample 
of BHB stars. The IXue et all (|2008f ) sample contains 
more than twice as many stars as our sample; however, 
our sample is deeper and contains twice as many stars 
with i? > 50 kpc. Neither sam ple fairly measur es the 
density profile of the halo. The IXue et al.l (|2008^ sam- 
ple is incomplete in color, magnitude depth, and spatial 
coverage. Our sample is complete in all dimensions, but 
we find an artificially shallow density profile - consis- 
tent with N{R) oc over the range 20 < i? < 70 
kpc (Figure [5]) - because we observe over a broader color 
range starting at go > 19 mag. A shallow density profile 
works to our advantage for the velocity dispersion analy- 
sis, however, because our stars sample greater distances 
at a greater relative density. 

3. HALO VELOCITY DISPERSION PROFILE 

Based on the distances and velocities of the sample 
stars, we calculate the velocity dispersion profile of the 
Milky Way halo. We use two independent methods to 
calculate velocity dispersion: a parametric method based 
on a Milky Way poten tial model, and a n on-pa rametric 
method based on the iDiaferio fc Gelled (|1997f) caustic 
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technique, originally developed to measure galaxy cluster 
mass profiles. Using these two methods provides a mea- 
sure of the systematic error introduced from clipping ve- 
locity outliers. We also discuss the systematic error from 
binary stars and disk stars in our dataset, and wc con- 
clude by comparing our velocity dispersion profiles with 
previous work. 

3.1. Computational Approach 

We use a Monte Carlo approach to model the distance 
to each star. We assume that photometric and velocity 
errors are Gaussian distributed, and that the underly - 
ing metallicity distribution is that of IXue et al.l (|2008D . 
We then derive the luminosity for each star by randomly 
drawing its color and metallicity from the observed distri- 
butions and comparing them to the main sequence and 
BHB tracks (Fig. [1]) described above. We sample the 
distributions 100 times per star. Using a Monte Carlo 
approach allows us to account for the non-Gaussian dis- 
tribution of luminosity estimates unique to each star. 

The resulting Monte Carlo catalog drawn from the 
observations produces the "cloud" of velocities and dis- 
tances shown in Figure [HI the distribution of points re- 
flects the uncertainties in the velocity measurements and 
the distance estimates of the stars. 

We calculate the velocity dispersion profile by group- 
ing stars into 5 or 6 unique bins in distance. The bins 
have sizes of 0.33i?, chosen so that the bins contain at 
least 70 stars. We require 50 stars to obtain a disper- 
sion with a formal statistical uncertainty of 10%. Bins 
with smaller occupation have dispersions systematically 
biased low resulting from small number statistics. Given 
the distance distribution of our sample, the occupation 
requirement constrains our velocity dispersion measure- 
ments to the region 15 < i? < 75 kpc. Stars with well- 
determined luminosities contribute their full weight to a 
single distance bin; stars with poorly-constrained lumi- 
nosities contribute less weight distributed over multiple 
bins. 

We use bootstrap re-sampling to calculate the uncer- 
tainty in the velocity dispersion measurement. Our pro- 
cedure is to draw random sets of 910 stars, with replace- 
ment, from the Monte Carlo catalog and re-calculate the 
velocity dispersion. We re-sample 10,000 times; the un- 
certainty is the standard deviation of the velocity disper- 
sions. 

3.2. Velocity Dispersion: Potential Model 

Our survey was designed to find unbound HVSs. Thus 
the first step in calculating the velocity dispersion profile 
is to clip the unbound stars from the sample. Defining an 
unbound star is difficult, however. The best estimate of 
the Galactic escape velocity is 550 ±50 km s~^ at the so- 
lar circle, base d on the highest v elocity stars in the solar 
neighborhood (jSmith et al.ll2007h . The escape velocity of 
the outer halo is more poorly constrained; extrapolating 
Vesc{R) to large R requires a pot e ntia l model. 

Here, we use the iKenvon et all (12008 ) potential model 
that is tied to observed Milky Way mass measurements. 
To establish VesciR) in this model, we drop a test par- 
ticle at rest from R = 500 kpc (i.e., halfway to M31) 
and calculate its tra jectory down to i? = kpc. The 
iKenvon et al.l (|2008D model predicts Vesc — 585 km 
at i? = 8 kpc in this definition of escape velocity, in good 




Fig. 6. — Upper panel: Monte Carlo distribution of distance 
and radial velocity for our sample. We represent each star by 100 
random draws from its parent distance and velocity error distribu- 
tion. We clip velocity outliers using two approaches: a parametric 
method based on a Milky Way potential model (dashed line), and 
a non-parametric method based on the .Diaferio fc Geller. ( 19971 ) 
caustic technique {solid line). Lower panel: binned velocity disper- 
sion profile from the two approaches. Errors are from bootstrap 
re-sampling. Linear least squares fits show an average —0.38 ±0.12 
km s~^ kpc~^ decline in velocity dispersion over 15 < i? < 75 kpc. 



agreement with the iSmith et al.l (|2007l ) observation. We 
fit a low-order polynomial to the calculated velocity as a 
function of distance and find 

VesciR) = -2.30xl0-^R^+0m88R'^-6.62R+Z km s-\ 

(3) 

where Z = 619 km s"\ valid for 15 < i? < 100 kpc. 

To avoid imposing an arbitrary mass on our velocity 
dispersion measurement, however, we re-normalize Equa- 
tion [3] to the observed envelope of negative velocity stars 
in our sample. Figure [S] shows the VesciR) relation that 
we use, with Z — 500 km s^^ (Eqn. [3]). Our assumption 
in changing Z to 500 km s^^ is that the velocity distribu- 
tion of halo stars extends up to the escape velocity, and 
that the most negative velocity stars - the stars falling 
in from the largest distances - provide the most robust 
measu re of escape velocity in our dataset. Percts et al] 
(|2009l) discuss this further in the context of HVSs. We 
test the effects of our choice of Vesc{R) by using a non- 
parametric approach to calculate velocity dispersion in 
the next section. In this section, we calculate the veloc- 
ity dispersion profile by using the shape of the VesciR) 
relation to provide a physically motivated means of clip- 
ping velocity outliers. 

We plot the resulting velocity dispersion of stars with 
velocities less than Vesc{R) for Z = 500 km s~^ in the 
lower panel of Figure \6\ Our Monte Carlo approach 
allows outliers that might otherwise be clipped to con- 
tribute a weight appropriate to their measurement un- 
certainties. A linear least squares fit to the velocity dis- 
persion profile finds a declining velocity dispersion. 



-0.30 ± 0.10)i? + (120 ± 4.6) km s" 



(4) 



valid over 15 < i? < 75 kpc. A higher-order fit does 
not significantly improve the residuals. The linear fit 
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has a standard deviation of 4.8 km s^ and a reduced 
of 0.7. A fixed velocity dispersion, by comparison, 

has a standard deviation of 7.6 km s~^ and a reduced 
of 1.5, poorer than the hnear fit but not statistically 

inconsistent. 

3.3. Velocity Dispersion: Caustic Method 

A non-parametric approach to determining the velocity 
dispersion profile leads to robust results based on fewer 
a-priori assumptions. Here we use the caustic technique. 
Because this technique is new to the stellar halo litera- 
ture, we begin with an introduct i on be fore discussing the 
mechanics and results. iDiafericI (|2009| ) provides a recent 
review of the technique. 

The caustic technique was originally developed to mca- 
sure galaxy clu ster mass profiles (Diafcrio & GcUcr 1997,; 
iDiaferid 1 1999( 1. Caustics are essentially escape velocity 
curves. The caustic technique is based on the distribu- 
tion of objects in phase space: the line-of-sight velocity 
V vs. projected distance r from the cluster center. In 
this plane, cluster members distribute in a characteristic 
trumpet shape with upper and lower borders indicating 
the escape velocity from the system. The amplitude A{r) 
of this distribution, which decreases with r, is a direct 
measure of the escape velocity from the system, inde- 
pendent of its dynamical state. The caustic technique 
only assumes that the spatial distribution of tracers is 
spherically symmetric. 

The caustic technique can be applied to any self- 
gravitating system. The technique is non-parametric and 
does not require that the tracers of the potential be in 
virial equilibrium. It is an effective tool for identifying 
bound members of the system defined by the upper and 
lower caustics (sharp declines in phase space density of 
tracers). For example. ISerra et al.l(|2009r ) apply the caus- 
tic technique to five dwarf spheroidals of the Milky Way. 

Applying the caustic technique to the halo star sam- 
ple requires a slight modification of the procedure used 
for galaxy clusters. Traditionally, the technique arranges 
sample objects in a binary tree according to their pair- 
wise projected binding energies (Diafcrio 1999). By walk- 
ing along the main branch of the tree, the technique de- 
termines the velocity dispersion which locates the caustic 
in the velocity diagram. This step does not apply to the 
Milky Way halo; the stars are not obviously clustered on 
the sky and the estimate of the pairwisc binding energy 
is inappropriate. Thus we apply the technique directly 
to the phase space diagram shown in Figure [6l 

The caustics are the curves satisfying the equation 
fq{r,v) = K, where fq{r,v) is the distribution in the 
phase space diagram, and k is the root of the equation 
{v^Bc)f^.R — 4cr^- The function 

{vlc).,R= r Al{rMr)dr/ ^{r)dr (5) 
Jo Jo 
is the mean caustic amplitude within R, (p{r) — 
J fq{'>',v)dv. R is not a free parameter in the stan- 
dard application of the technique, but here we chose 
R — 30 kpc. Our results are totally insensitive to this 
parameter: varying R in the range 20-100 kpc varies the 
final number of the halo members by at most five and 
does not change the velocity dispersion profile. 

We use an iterative approach to estimate the velocity 
dispersion profile. First we compute the velocity dis- 



persion of the total sample and locate the caustics that 
enable some interloper removal. We then compute a new 
velocity dispersion with the member stars, locate new 
caustics and remove further interlopers. We proceed un- 
til no star is identified as an interloper. This procedure 
requires only four steps and on average yields 838 (of 
910) final halo members and an overall velocity disper- 
sion of 99 km s~^. The final set of caustics are drawn 
with the solid line in the upper panel of Figure [H 

The open triangles in the lower panel of Figure [S] show 
the caustic velocity dispersion profile. It is visually ap- 
parent that the caustic technique measures the velocity 
dispersion from the well-sampled core of the velocity dis- 
tribution. This approach results in a reliable velocity 
dispersion profile, but yields a velocity dispersion that is 
systematically 10% smaller than found by our paramet- 
ric Vesc{R) method. We conclude that clipping velocity 
outliers with the caustic technique changes the observed 
amplitude, but not the observed slope, of the velocity 
dispersion profile. 

A linear least squares fit to the caustic velocity disper- 
sion profile finds a declining velocity dispersion, 

ay = (-0.45 ± 0.16)i? + (110 ± 6.0) km (6) 

valid over 16 < i? < 64 kpc. The linear fit has a standard 
deviation of 4.8 km s~^ and a reduced of 1.0. A fixed 
velocity dispersion, by comparison, has a standard devi- 
ation of 7.9 km s~^ and a reduced of 2.8. In this case, 
a constant velocity dispersion is a significantly poorer fit 
to the data. 

3.4. Possible Systematics 

There are at least two contaminants that may sys- 
tematically affect our velocity dispersion profile: binary 
stars, which increase the velocity dispersion, and disk 
stars (i.e. white dwarfs), which decrease the velocity dis- 
persion. We investigate how these contaminants may 
alter the observed velocity dispersion profile. 

Binary stars are unlikely to change the velocity disper- 
sion profile. Both BHB stars and 1.5 M0 blue stragglers 
have stellar radii around 2.5 Kq. Equal mass pairs of 
such stars must have semimajor axes of at least 6.5 R0 to 
avoid Roche lobe overflow. Thus the most compact pos- 
sible binary system has velocity semi- amplitude of 100 
km s^^. Assuming that half of the targets are binaries 
with a log-normal distribution of semimajor axes, we ex- 
pect 30 binaries in our dataset with wsini >50 km s~^. 
This estimate is generous given that BHB stars have re- 
cently evolved through the red giant phase and are thus 
unlikely to have close companions. In any case, binaries 
will be observed at a random orbital phase and binned 
with 70-200 other stars to compute a velocity disper- 
sion. We propagate our simulated distribution of v sinj's 
through our Monte Carlo catalog and find a negligible 
change in the velocity dispersion profile; the uncertainty 
in the velocity dispersion slope remains 30%. 

Disk stars pose a greater threat to the velocity dis- 
persion profile. Disk stars have a systematically lower 
velocity dispersion than halo stars and may also have a 
non-zero mean velocity because of the longitude depen- 
dence of the Sun's circular velocity correction. The SDSS 
imaging survey from which we draw our candidates does 
not uniformly survey the sky across all longitudes. 
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Our greatest concern is white dwarfs, which are intrin- 
sically faint and overwhelmingly appear at faint magni- 
tudes in our survey. Inserting white dwarfs into our halo 
sample thus reduces the velocity dispersion observed at 
large distances. Observationally, 15% of our survey tar- 
gets are white dwarfs, mostly found at low latitudes. If 
we insert a 1% white dwarf contamination into our Monte 
Carlo catalog, the resulting velocity dispersion profile 
steepens by 10%. Fortunately, white dwarf s are read- 
ily ide ntified by their broad Balmer lines (see iKilic et all 
l2007al) . Visual inspection of our spectra reveals no white 
dwarf contaminants. 

Binaries also have observational constraints. We ob- 
tained repeat observations for every star with \vrf\ > 
300 km s~^ and found only one star, a star near —300 
km s~^, to exhibit radial velocity variation. We thus ex- 
cluded this star from the halo sample. We conclude that 
binaries and white dwarfs are unlikely to significantly in- 
fluence the velocity dispersion profile measured from our 
halo star sample. 

3.5. Comparison with Previous Results 

There are few measurements of the Milky Way veloc- 
ity dispersion profile at lar ge R. The most comparable 
measu r ements are [Batt a glia et al.l ([2005) and Xue ct al. 
([20081 ). iBattagha etUl (2005') measure the velocity dis- 
persion profile based on a sample of 9 satellite galaxies, 
44 globular clusters, 58 red giants, and 130BHB stars 
that span 10 < < 140 kpc. lXue et al] ([200g ) measure 
the velocity dispersion profile based on a sample of 2401 
BHB stars that span 5 < i? < 60 kpc. With the excep- 
tion of a handful of shared BHB stars, the samples are 
independent of one another. 

Figure [7| compares the velocity dispersion profile and 
radial number d istribu tion of our sam p le wit h those of 
IBattagha eFall (|20050 and 'Xu e et"al] ((2008). Beyond 
R > 50 kpc, our sample con tains a factor of 8 more stars 
than [Battaglia et al.l ([2005| ) and a factor of 2 more stars 



than [Xue et al.[ ([2008D . Remarkably, the velocity disper- 
sions measured by the three samples are consistent at the 
1.5-cr level. The three samples also observe a statistically 
similar decline in velocity dispersion with distance. The 
major difference between the samples is the significance 
of the observed decline in velocity dispersion with dis- 
tance. 

A linear least squares fit to the [Xue et all ([2008D data 
gives a —0.15 ± 0.14 km s~^ kpc~^ decline in velocity 
dispersion over 10 < i? < 60 kpc. The linear fit has a 
standard deviation of 6.3 km s~^ and a reduced of 
1.5. A fixed velocity dispersion has a very similar stan- 
dar d deviation of 6.4 km s~^ and reduced of 1.9. Thus 
the [Xue et al.[ ([2008D measurements cannot formally dis- 
criminate between a constant velocity dispersion and a 
declini ng velocity disp e rsion. 

The [Battaglia et al.l ([2005[ ) data, on the other hand, 
exhibit a steeper —0.58 ± 0.11 km kpc""'^ decline in 
velocity dispersion over 10 < i? < 100 kpc. We exclude 
their last bin because it contains only 3 objects. The 
linear fit has a standard deviation of 8.6 km s~^ and a 
reduced of 0.85. 

The weighted average of all the samples yields a 
—0.38 ± 0.12 km s~^ kpc~^ decline in velocity disper- 
sion over 10 < i? < 100 kpc. We obtain the same result 
if we average only our own two velocity dispersion mea- 
surements. The Milky Way velocity dispersion profile is 
not linear, of course, but three independent sets of ob- 
servations are in statistical agreement: the Milky Way 
radial velocity dispersion drops from a ~ 110 km at 
R=lb kpc to CT ~ 85 km s'^ at i? = 80 kpc. 

4. CONCLUSIONS 

The mass and mass distribution of the Milky Way are 
fundamental parameters because they link directly to 
theoretical models. We use a spectroscopic sample of 
910 halo stars derived from our HVS survey to measure 
the velocity dispersion profile of the Milky Way. The 
stars are 74% BHB stars and 26% blue stragglers. We 
estimate luminosities using stellar evolution tracks for 
metal poor main sequence stars and BHB stars. Because 
of the non-Gaussian distribution of luminosity estimates, 
we use a Monte Carlo approach to calculate velocity dis- 
persion and its uncertainty. 

We calculate the velocity dispersion profile in two ways: 
a parametric method based on a VesciR) model, and a 
non-parametric method based on the caustic technique 
originally developed to measure galaxy cluster mass pro- 
files. Comparing the two methods provides a measure 
of the systematic uncertainty arising from the clipping 
of outliers in velocity. The velocity dispersion from the 
caustic method is 10% smaller than the Vesc{R) method, 
but both methods identify a similar decline in velocity 
dispersion with distance: —0.38 ± 0.12 km s^^ kpc^^, 
valid for 15 < i? < 75 kpc. 

Our sample c onta ins a factor of 8 more stars than 
Battaglia et al.l ([200 5) and a factor of 2 more stars than 
Xue et air i 20081 ) at i? > 50 kpc. The velocity dispersion 
profiles observed by these independent datasets are con- 
sistent at the 1.5-CT level, and have an average velocity 
dispersion slope identical to our result. Remarkably, no 
matter what tracers are used, observers find the same 
halo velocity dispersion profile. 

The velocity dispersion profile is a basis for measuring 
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TABLE Al 
DATA TABLE 



RA Dec go (Tg {u-g)o ^{u-g) id - 1-}o ^{g^r) ""helio I b V^f Mg^BHB CBH 

hrs deg mag mag mag mag mag mag km km deg deg km mag ma 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15 



0:02:05.713 31:18:50.23 20.434 0.026 0.778 0.095 -0.277 0.040 -209.3 14.0 110.717 -30.385 -34.3 1.17 0. 

0:02:33.817 -9:57:06.85 18.434 0.021 0.753 0.040 -0.328 0.040 -87.7 9.7 86.763 -69.316 -14.8 1.31 0. 

0:04:36.491 -9:57:19.48 19.834 0.018 0.922 0.103 -0.167 0.032 -172.6 35.0 87.977 -69.584 -100.7 0.69 0. 

0:05:28.141 -11:00:10.07 19.141 0.042 1.007 0.081 -0.275 0.047 -115.9 10.9 86.890 -70.586 -47.8 0.94 0. 

0:07:52.013 -9:19:54.32 17.302 0.017 1.016 0.036 -0.276 0.039 -114.5 9.9 90.855 -69.443 -42.2 0.91 0. 

0:12:26.890 -10:47:54.56 18.898 0.025 1.006 0.064 -0.321 0.038 -128.2 10.5 91.741 -71.269 -62.8 0.94 0. 

0:23:53.294 -1:04:46.40 18.199 0.016 0.749 0.042 -0.255 0.025 19.7 9.8 107.552 -63.125 109.0 1.15 0. 

0:29:31.158 15:39:40.20 19.069 0.024 1.057 0.069 -0.270 0.038 22.3 35.0 115.201 -46.881 153.4 0.88 0. 

0:36:40.570 -11:11:25.02 17.421 0.018 0.778 0.028 -0.304 0.031 32.5 9.8 109.940 -73.689 84.1 1.26 0. 

0:39:06.749 24:09:05.62 19.348 0.018 0.948 0.081 -0.255 0.030 -130.0 35.0 119.332 -38.634 15.0 0.94 0. 



Note. — Tablc fATl is presented in its entirety in the eleetronic edition of the Astrophysieal JournaL A portion is shown here for guidance and content. 
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For further progress in measuring the Milky Way ve- proved this paper. This project makes use of data prod- 
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at distances i? > 50 kpc. It is difficult to find i? > 50 by the Astrophysieal Research Consortium for the Partic- 
kpc tracers because of the steep decline in the density of ipating Institutions. This research makes use of NASA's 
the stellar halo. It is also very difficult for proper mo- Astrophysics Data System Bibliographic Services. AD 
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APPENDIX 
DATA TABLE 

Table lAll presents the 910 stars used here. We provide the observed positions, magnitudes, and velocities plus our 
derived luminosities, distances, and BHB likelihood. The table columns are: (1) RA (J2000), (2) Dec (J2000), (3) 
dereddened SDSS .go magnitude, (4) magnitude error, (5) dereddened SDSS {u — g)o color, (6) (u — g)o error, (7) 
dereddened SDSS (g — r)o color, (8) {g — r)o error, (9) hehocentric radial velocity Vheiio, (10) velocity error, (11) 
Galactic longitude I, (12) Galactic latitude b, (13) Galactic rest f rame velocity Vrf, as defined in Equation [1] (14) 
BHB absolute magnitude Mg^sHB derived from lDotter et al.l ([IDOS") for [Fe/H]=-1.9, (15) magnitude error, (16) BHB 
Galactocentric distance Rbhb, (17) blue straggler absolute magnitude Mg,BS derived from Girardi et al. (2004) for 
[Fe/H]=— 1.9, (18) magnitude error, (19) blue straggler Galactocentr ic d istance Res, and (20) the star's likelihood of 
being BHB, < Jehb < 1, based on stellar colors and spectra. Table lAll is available in its entirety in machine-readable 
form in the online journal. A portion of the table is shown here for guidance regarding its form and content. 
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